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ABSTRACT 

^  The  determination  of  the  distribution  of  airborne  toxic  particles  as  a  function 
of  the  aerodynamic  diameter  provides  important  information  as  weil  as  criteria  for  the 
definition  of  hazard  as  applied  to  levels  of  airborne  contamination.  This  is  because  the 
aerodynamic  particle  size  distribution  embodies  the  information  related  to  particle  den¬ 
sity,  diameter,  shape  factor  and  slip  correction  that  is  critical  for  the  characterization  of 
particle  motion  in  settling  and  impaction  and  it  is  these  motions  that  are  responsible 
for  particle  deposition  in  the  respiratory  tract  and  particle  collection  in  aerosol  sampling 
devices.  For  a  given  definition  of  hazard  based  on  some  parameter  related  to  the 
aerodynamic  size  distribution,  this  paper  develops  a  statistical  sampling  error  model  for 
the  parameter  that  is  based  on  the  Poisson  process.  Given  that  an  appropriate  sampling 
program  has  been  designed  for  the  measurement  of  the  size  distribution-related  parameter 
with  the  aerodynamic  particle  size  analyzer,  this  paper  proceeds  to  the  derivation  of  an 
optimum  detection  algorithm  for  the  detection  of  a  signal  aerosol  sequence  in  a  set  of 
J  aerosol  samples  with  a  common  background.  The  detection  algorithm  is  based  on  the 
generalized  likelihood  ratio  test  in  which  the  received  count  associated  with  the  aerosol 
sample  is  modeled  as  a  Poisson  distributed  random  variabla^rformance  analyses  of 
the  resulting  algorithm,  based  on  the  probability  of  detection  (Pp)\ersus  signal-to-noise 

ratio  for  several  given  fixed  false  alarm  probabilities  (Pp^),  are  pres^ted. 
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RESUME 


La  determination  de  la  distribution  des  particules  en  suspension  dans  I’air 
en  fonction  de  leurdiametre  aerodynamique  fournit  d’importantes  informations  ainsi  que 
des  criteres  qui  permettent  de  definir  les  risques  correspondent  aux  degr^s  de  contamina¬ 
tion  du  milieu  atmospherique.  En  effet,  la  distribution  de  lataille  aerodynamique  renferme 
des  informations  sur  la  masse  volumique,  le  diametre,  le  facteur  de  forme  et  la  correc¬ 
tion  de  glissement  des  particules,  qui  sont  essentielles  pourcaracteriser  le  mouvement 
des  particules  au  cours  de  leur  depdt  et  de  leur  impact;  ce  mouvement  est  ^galement 
responsable  du  d^pdt  des  particules  dans  les  voies  respiratoires  et  de  leur  pi^geage  dans 
les  dispositifs  de  pr^l^vement  d’aSrosols.  Pour  une  definition  donn^e  du  risque  deter¬ 
mine  e  partir  d’un  parametre  quelconque  lie  e  la  distribution  de  la  taille  aerodynamique, 
on  eiabore  dans  cette  communication  un  moddle  statistique  base  sur  la  loi  de  Poisson, 
permettant  de  determiner  I'erreur  d’echantillonnage  pour  le  parametre.  A  I’aide  d’un  pro¬ 
gramme  de  preievement  approprie  congu  pour  mesurer,  avec  I’analyseur  de  tallies 
aerodynamiques,  le  parametre  lie  e  distribution  de  la  taille,  on  deduit,  dans  cette  com¬ 
munication,  I’algorithme  optimal  permettant  de  detecter  une  sequence  d’aerosols  dans 
un  ensemble  de  jechantillons  d'aerosols  qui  presentent  les  memes  caracteristiques  de 
base.  L’algorithme  de  detection  est  base  sur  un  test  generalise  du  rapport  des 
vraisemblances,  dans  lequel  la  valeur  regue  associee  k  rechantillon  d’aerosol  est 
modeiisee  comme  une  variable  aieatoire  qui  obeit  a  la  loi  de  Poisson.  On  presente  les 
resultats  des  analyses  de  performance  de  I’algorithme  ainsi  obtenu,  bases  sur  la  pro- 
babilite  de  detection  (P^)  en  fonction  du  rapport  signal/bruit  pour  plusieurs  valeurs  don- 
nees  de  probabilites  de  fausse  alarme  (Pp^)- 
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I.  INTRODUCTION 

The  aerodynamic  particle  sizer  (APS),  Model  3300  (TSI  Incorporated)  can  provide  a  real¬ 
time  measurement  of  the  particle  size  distribution  according  to  the  aerodynamic  diameter  [1].  This 
sampling  instrument  determines  the  aerodynamic  diameter  of  the  individual  aerosol  particles  in 
the  sample  by  measuring  the  transit  time  of  the  peirticles  between  two  spots  generated  by  a  laser 
velocimeter  that  employs  a  2  mW  polarized  He-Ne  laser  as  the  light  source. 

The  APS  consists  of  two  basic  modules:  (1)  the  sensor  module  and  (2)  the  signal  processing 
and  microcomputer  module.  The  main  parts  of  the  sensor  module  are  the  accelerating  nozzle  and 
the  two  spot  laser  velocimeter.  This  module  brings  the  aerosol  sample  into  an  outer  accelerating 
orifice  and  focusses  the  individual  particles  through  a  dual  beam  laser  formed  by  splitting  a  focussed 
laser  beam,  on  the  basis  of  polarization,  using  a  calcite  plate.  The  beams  are  then  focussed  using 
a  cylindrical  lens  to  produce  two  flat  beams  of  rectangular  cross-section  just  downstream  of  the 
nozzle  orifice.  As  the  aerosol  particle  passes  through  these  two  beams,  it  triggers  a  pair  of  electrical 
pulses  whose  temporal  separation  is  accurately  measured  using  a  high  speed  digital  clock.  The 
conversion  of  this  information  into  an  aerodynamic  particle  size  is  directed  by  the  signal  pro  ^ing 
module.  A  multichannel  accumulator  (MCA)  is  used  to  record  the  transit  times  of  all  aerosol 
particles  and,  at  the  end  of  the  prescribed  sampling  period,  a  microcomputer  reads  each  channel  of 
the  MCA,  translates  the  channel  numbers  to  aerodynamic  particle  sizes,  and  displays  the  information 
as  a  discrete  spectrum  or  histogram  consisting  of  48  size  intervals  (i.e.,  bins)  spanning  the  aerosol 
diameter  range  from  0.5-15  /tm. 

Currently,  there  is  active  interest  in  the  development  of  signal  proce.ssing  algorithms  for  the 
deUrction  of  weak  aerosol  target  signals  in  a  heavy  clutter  aerosol  background  using  an  APS  as  the 
sensor  detection  system.  A  theoretical  investigation  of  the  performance  characteristics  of  the  APS, 
when  utilized  in  the  detection  mode,  is  undertaken  in  this  paper.  Towards  this  objective,  a  stochastic 
sampling  error  model  and  a  signal  detection  algorithm  based  on  the  application  of  the  classical 
generalized  maximum  likelihood  ratio  test,  are  developed  for  the  APS.  This  paper  is  organized  as 
follows.  In  Section  II,  a  statistical  sampling  error  model  for  the  APS,  based  on  the  Poisson  process, 
is  developed  and  analyzed.  A  signal  detection  algorithm  which  applies  the  generalized  mEtximum 
likelihood  ratio  test  to  the  Poisson  process  is  developed  in  Section  III.  In  order  to  analyze  the 
performance  of  the  detection  algorithm,  probabilities  of  false  alarm  and  detection  as  a  function  of 
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the  aerosol  target-to-clutter  signal-to-noise  ratio  (SNR)  are  also  presented  in  Section  III.  Section  IV 
provides  an  application  of  the  statistical  sampling  error  model  and  the  signal  detection  algorithm  to 
two  explicit  numerical  examples  involving  the  detection  of  a  transient  aerosol  target  signal  with  an 
unknown  arrival  time  in  a  background  clutter  aerosol  environment.  Finally,  conclusions  are  listed  in 
Section  V. 

II.  A  STATISTICAL  SAMPLING  ERROR  MODEL  FOR  THE  APS 

The  measurement  of  the  aerodynamic  particle  size  distribution  function  (PSDF)  or,  for  that 
matter,  any  functional  of  the  PSDF  is  necessarily  subject  to  an  intrinsic  variability  that  is  due  both  to 
statistical  sampling  errors  and  to  natural  physical  variations  which  are  not  readily  separable  from  the 
statistical  errors.  For  a  single  sampling  instrument  such  as  the  APS,  the  statistical  sampling  errors 
can  be  reduced  by  the  design  of  an  appropriate  sampling  method.  This  section  is  concerned  with 
the  estimation  of  the  sampling  errors  associated  with  the  operation  of  the  APS  with  the  objective  of 
specifying  the  required  sample  size  and  observation  times  that  must  be  adopted  in  order  to  bound 
such  errors. 

The  sampling  error  model  for  the  APS  Is  based  on  the  Poisson  model  which  provides  the 
information  about  the  statistical  properties  of  the  random  point  process  associated  with  the  fluc¬ 
tuations  of  any  size  distribution-related  variable  as  a  consequence  of  the  stochastic  fluctuations  of 
aerosol  particle  numbers  in  a  given  size  range  in  the  sampling  volume.  It  is  asserted  that  the  particle 
number  fluctuations  in  a  given  aerodynamic  diameter  range  can  be  adequately  modeled  by  a  Poisson 
random  variable  [2]  with  the  probability  A:  =  0, 1,2, . . .  particles  occurring  given  by 

A* 

Pit  =  —  exp(-A),  A>0, 

where  A  is  the  Poisson  parameter. 

With  the  assumption  that  the  Poisson  distribution  properly  describes  the  random  fluctu¬ 
ations  of  aerosol  particle  numbers  of  a  given  aerodynamic  diameter  range  in  the  sampling  volume 
Vi,  the  variability  of  any  size  distribution-related  variable  can  be  derived  as  follows.  For  any  given 
function  f(D)  of  the  aerodynamic  diameter  D,  any  particular  measurement  of  the  observed  value 
of  the  function  due  to  aerosol  particles  in  the  size  interval  D  to  D  +  dD  in  the  sampled  volume  V„ 
is  given  by  dF  =  f{D)n{D)  V,dD,  where  n{D)  is  a  particular  realization  of  the  aerodynamic  par¬ 
ticle  size  distribution.  Consequently,  it  can  be  concluded  that  the  ensemble  mean  of  the  stochastic 
functional  F,  due  to  the  contribution  of  all  aerosol  particles  within  the  aerodynamic  diameter  range 
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[Di,  £>2]  in  the  sampling  volume  V„  is  given  by 

fDj 


{F)=  f  f(D){niD))V,dD 
JDt 

=  f  'f(D)N(D)V,dD, 

JDi 


(1) 


where  the  ensemble  average  of  any  quantity  is  indicated  by  angled  brackets,  (  ).  It  is  worthwhile 
pointing  out  that  the  ensemble  mean  number  of  aerosol  particles  with  dieunetcrs  D  to  D  +  dD  in 
volume  \i  is  given  by  N(D)  V,dD  =  {n(D))  V,dD.  In  other  words,  the  aerosol  particles  sampled  in 
the  volume  If  are  Poisson  distributed  about  the  ensemble  mean  particle  size  distribution,  N{D). 
Similarly,  the  variance  of  F  can  be  written  as 

fOi 

\ht(F)  =  aj.=  /  /^(D)  Var{n(D))  V,  dD. 

Jd, 

Now,  recall  that  the  mean  and  variance  of  a  Po'sson  distribution  are  equal.  Since  it  has  been  assumed 
that  n(D)  at  any  fixed  value  of  D  is  a  Poisson  random  variable,  it  follows  that  Var(n(D))  =  o^d)  — 
{n(D))  =  N(D)  so 

^ f\D)N{D)V,dD.  (2) 

JDt 


In  light  of  Eqs.  (1)  and  (2),  the  relative  fluctuation  intensity  of  F  defined  as  Ip  =  ap/iF), 
can  be  expressed  as  follows: 

ro,  \  >/2 

(  / 

ap 


Ip  = 


(^l\\D)N{D)V,dD^ 

f{D)N{D)V,dD 

Jd, 


(3) 


Since  several  particulate  proportits  are  proportional  to  the  diameter  of  the  particles  raised  to  an 
appropriate  power,  it  is  useful  to  coasider  functions  /(D)  with  the  form  f{D)  =  c„D",  whore  c„  is 
a  constant  shape  factor  whose  precise  form  is  determined  only  by  the  specific  value  of  n.  Examples 
of  such  functions  include  the  following:  J{D)  =  1  (i.c.,  n  =  0  and  Q)  =  1)  for  which  the  associated 
property  is  the  expected  value  (ensemble  average)  of  the  total  number  of  aerosol  particles  within 
the  size  interval  {Di,  £>2)  in  the  sampled  volume;  /(D)  =  ttD*  (i.e.,  n  =  2  and  C2  =  tt)  for  which  the 
associated  property  is  the  equivalent  surface  area  of  the  aerosol  particles;  and,  /(D)  =  7rpD’/6  (i.e., 
n  =  3  and  C3  =  7rp/6  where  p  is  the  particle  density)  for  which  the  associated  property  is  the  mass 
of  the  aerosol  particles.  Insertion  of  this  specialized  form  for  /(D)  into  Eq.  (3)  yields 

aDj  \  1/2 

D'^’'N{D)V,dDj 


T  =  - 

~  (^n> 


fDt 

/  D"N(D) 
Jd, 


(4) 


V,dD 
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Since  the  with  moment  of  the  particle  diameter  about  the  origin  D  =  0  over  the  size  interval  [Dj,  Ih] 
is  defined  as 

fOi 

/  D'"N(D)  dD 

^  - < 

/  N{D)  dD 

JDi 

it  follows  easily  that  the  relative  fluctuation  intensity  (cf.  Eq.  (4))  can  be  written  in  the  following 
form; 

-  (F„)  ^  ^  (D")  ’ 

where  (JV),  the  expected  total  number  of  particles  in  the  size  diameter  range  [Di,  Dj]  in  the  sampled 
volume  Vf,  is  given  by 

(N)  =  f  N{D)  Vi  dD. 

JDi 

It  is  informative  to  write  Eq.  (5)  as 

^  =  N,  'l\  (flu) 

\^n) 

where  Ne  is  an  effective  aerosol  particle  number  which  is  defined  by 

N,  =  k(N)  (66) 


In  light  of  this  development,  it  is  seen  that  the  relative  fluctuation  intensity  of  any  particle  size 
distribution- related  variable  F„  is  proportional  to  the  inverse  square  root  of  the  effective  total  number 
of  aerosol  particles,  Ng.  The  effective  total  number  of  particles,  is  directly  proportional  to  the 
expected  value  of  the  total  number  of  aerosol  particles,  {N),  with  the  constant  of  proportionality, 
K,  involving  the  moments  of  the  particle  size  distribution  as  per  Eq.  (6c). 

In  order  to  provide  some  concrete  results  concerning  the  statistical  sampling  error  model 
derived  above,  it  is  useful  at  this  point  to  choose  an  explicit  form  for  the  aerodynamic  particle  size 
distribution  function,  N{D).  For  this  purpose,  it  is  assumed  that  the  aerosol  particle  size  distribution 
can  be  adequately  approximated  by  a  three  parameter  gamma  size  distribution  represented  by 


where  a>0,  /?>0,  D>0  and 


fOO 

r{a)  =  I  t“~*cxp(t)dt 

Jo 
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is  the  gamma  function.  Here,  a  is  the  shape  parameter,  ^  the  scale  parameter,  and  N'  the  particle 
number  density.  It  is  noted  that  the  gamma  distribution  provides  an  attractive  compromise  between 
the  widely  used  log-normal  and  exponential  distributions  that  has  been  traditionally  assumed  for 
the  aerosol  size  distribution.  Indeed,  measurements  of  Berger,  Melice  and  Demuth  [3]  to  test  the 
applicability  of  the  log-normal  and  gamma  distribution  on  an  extensive  data  set  for  suspended 
particulates  in  Gent,  Belgium  have  demonstrated  conclusively  that  the  gamma  distribution  was 
superior  to  the  log-normal  distribution  for  modeling  the  suspended  particulate  levels.  For  the  gamma 
size  distribution,  N{D),  the  relative  fluctuation  inteiKsity  />;  expressed  in  ICq.  (4)  reduces  to  the 
following  form: 

1  1  [Pj^n  4-  g,  Ih/fi)  -  F(2«  -h  a,  Dj/fi)] 

{Fn)  P{n+a,D2l(i)-P{n+a,Dxll3)  ’ 

where  P(a,  x)  is  the  incomplete  gamma  function  defined  as 

Similarly,  it  can  be  readily  shown  that  the  parameter  k  for  the  gamma  size  distribution  can  be 
represented  by 

^ _ [P{n+a,D2/(i)  -  P{n+  a.Pi//?)]^ _ 

[P{a,  D2I13)  -  P(a,  Dxll3)]  [P(2n-f  a,  -  P(2n-f-  cr,  Dx/P)] '  ^  ^ 

For  a  fixed  size  interval  [D\ ,  Dt],  the  pwameter  k  is  a  function  only  of  jt,  the  exponent  related 
t  o  the  size  distribution-related  variable  F„  atid  of  a  and  /#,  the  gamma  shape  and  scale  parameters, 
respectively.  A  plot  of  k  =  1/k  =  {N)/Nc  as  a  function  of  7i  for  a  =  1,2  and  0  =  1.0  fim,  Dx  =  0.5 
/tm  and  D2  =  15.0  /tm  is  shown  in  Figure  1.  It  is  evident  from  this  plot  that  as  the  exponent 
n  increases,  thereby  providing  increasing  weight  to  the  aerosol  particles  with  the  larger  diameters, 
K  (and,  therefore  the  effective  particle  number,  Nc,  as  per  Eq.  (6b))  decreases  very  rapidly,  thus 
increasing  the  relative  fluctuation  intensity  of  F„.  Intuitively,  this  is  due  to  the  smaller  number  of 
aerosol  particles  which  contribute  to  the  major  fraction  of  F„  as  n  increases.  Furthermore,  note  that 
K  =  1  when  n  =  0,  with  the  result  that  in  this  case,  Eq.  (6a)  reduces  to  on/{N)  =  The 

latter  relationship  is  a  well-known  result  for  the  relative  fluctuation  intensity  of  the  total  number  of 
aerosol  particles,  N,  in  the  sampled  volume 

By  application  of  Fzp  (7),  it  is  possible  to  construct  curves  for  the  normalized  relative 
fluctuation  intensity  for  various  size  distribution-related  variabhs.  Towanls  t.his  end,  we  consider 
the  size  interval  (/)i,/^2]  =  [0.5,15.0]  /irn  which  is  the  operating  particle  size  range  for  the  APS. 
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With  this  operating  constraint,  curves  for  the  normalized  relative  fluctuation  intensity,  deflned  as 
(Af*  have  been  calculated  as  a  function  of  the  gamma  size  distribution  shape  parameter  a 

and  are  displayed  in  Figures  2,  3,  and  4  for  /S  =  1.0  /tm  and  n  =  0,2,3,  respectively.  These  curves 
distinctly  show  the  variation  of  the  normalized  relative  fluctuation  intensity  of  F„  with  the  shape 
of  the  underlying  number  size  distribution  function  as  embodied  in  the  gamma  shape  parameter  a. 
To  provide  further  illustration  of  the  effects  of  the  number  size  distribution  shape  on  the  relative 
fluctuation  intensity,  we  consider  the  ratio,  R,  deflned  as  the  ratio  of  If„  using  the  gamma  size  dis¬ 
tribution  to  that  using  the  exponential  size  distribution  (to  which  the  gamma  distribution  converges 
as  Q!  — ►  1).  In  light  of  Eq.  (7),  the  ratio  R  can  be  expressed  as 

^  _  Jf„  (gamma) 

(exponential) 

^  [F(2n  -f  o.  Eh!  13)  -  P(2n  -i-  a,  Di/fi)]  [P(n  -h  1 ,  D;//?)  -  P(n  -H,  Di/H)] 

[F(2n-h  1,D2/^)  -  P(2n-(- 1,  D, //#)]*/' [P(n-f  a,  0.113)  -  F(n a,  Dj//?)]  ‘ 

Plots  of  the  ratio  R  with  changing  a  for  the  cases  of  n  e<)ual  to  0,  2  and  3  are  display«l  in  Figure  5 
for  (3  =  1.0  /tm  and  in  Figure  6  for  /I  =  3.0  /trn.  Tho.se  plots  provide  a  visual  intlicalion  of  the 
distribution  shape  effects  on  the  relative  fluctuation  intensity  of  the  sampled  variable  Fn,  basotl 
on  differences  between  the  two  sampling  error  models.  The  most  important  point  to  note  in  these 
plots  is  that  there  can  be  a  change  in  the  relative  fluctuation  Intensity  of  the  variable  F„  that  arises 
simply  because  of  a  variation  in  the  size  distribution  shap>e.  For  example,  when  a  gamma  size 
distribution  having  an  actual  shape  parameter  a  =  0.1  {fi  —  1.0  /tm)  is  sampled  with  an  APS  and 
the  relative  fluctuation  intensity  of  the  number  concentration  (i.e.,  n  =  0)  in  the  diameter  range 
\Di,D^  =  [0.5,15.0]  /im  is  subsequently  computed  using  an  exponential  model  (i.e.,  a  — »  1),  then 
the  actual  relative  fluctuation  intensity  of  the  number  concentration  in  this  case  is  about  5  times  the 
one  computed.  This  implies  that  there  is  an  underestimation  by  about  5  times  the  actual  relative 
fluctuation  intensity  value,  which  can  arise  simply  by  neglecting  the  effects  of  the  distribution  shape. 

As  a  final  point,  it  should  be  noted  that  the  statistical  sampling  error  motlel  dovelofxxl  in 
this  section  can  be  used  to  compute  the  sampling  errors  for  a  large  number  of  variables  relaUxI  I.o 
the  aerodynamic  particle  size  distribution.  In  particular,  the  model  can  be  applied  to  determine  the 
minimum  sampling  duration,  that  would  be  required  for  operation  of  the  APS  at  a  given  aerosol 
sampling  flowrate,  Q„  in  order  that  the  sampling  error  of  a  given  variable  of  interest  (e.g.,  F„)  be 
smaller  than  or  equal  to  some  prescribed  Jf„-  As  a  simple  example,  to  achieve  a  sampling  error  of 
one  percent  for  the  total  aerosol  particle  number  using  the  APS  with  Dj  =  0.5  /im  and  Z>j  =  15.0 
fim  for  an  underlying  size  distribution  function  with  a  =  \  .b,  0  =  2.5  /tm,  and  N*  =  6.0  cm"*  will 
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require  a  minimum  sampling  volume  of  =  1800  cm®  (cf.  Eq.  (7)).  For  the  Al’S  operating  at  an 
aerosol  flowrate  of  Qt  =  83.3  cm®  s~\  this  requires  a  sampling  duration  of  ^  =  22  s. 

III.  DETECTION  ALGORITHM  AND  PERFORMANCE  ANALYSIS 

In  this  section,  a  detector  based  on  the  parameter  F„,  whose  statistical  sampling  errors  were 
analyzed  in  the  previous  section,  is  proposed.  Towards  this  objective,  consider  the  APS  sampling 
for  duration  ^  at  an  aerosol  flowrate  Q,  to  produce  count  data  as  the  set  of  numbers  {F„(j),j  = 
1, 2, . . . ,  J},  viz.,  these  numbers  are  the  observed  values  of  the  parameter  F„  due  to  aerosol  particles 
in  the  size  interval  |Di,D2]  in  the  sampled  volume  V,  =  Q,tg.  These  counts  are  random  variables 
and,  as  in  the  previous  section,  it  is  assumed  that  the  count  F„(j)  during  the  jth  sampling  period 
(reception  or  measurement)  is  a  Poisson  random  variable  that  assumes  the  values  A:  =  0, 1, 2, . . .  with 
probability  Pt  =  (A*exp(— A))//:!  (A  >  0).  The  counts  obtained  during  disjoint  sampling  periods 
(i.e.,  exposures)  are  assumed  to  be  independent  random  variables. 

For  the  purposes  of  detection,  it  is  assumed  that  the  aerosol  sample  is  the  sum  of  a  target 
“signal”  aerosol  whose  presence  is  to  be  detected  and  a  background  “clutter”  aerosol  whose  presence 
is  unavoidably  received  from  the  natural  aerosol  component  present  in  the  atmosphere.  In  this  case, 
each  count,  F„(j),  will  be  a  Poisson  random  variable  with  parameter  A^*^0)  =  ^0)  +  where 

refers  to  the  situation  when  the  signal  aerosol  to  be  detected  is  actually  present  in  the  sample, 
a  situation  that  is  referred  to  as  hypothesis  Hi.  When  the  signal  aerosol  is  absent  in  the  sample — a 
situation  referred  to  as  hypothesis  Hq — the  count  data  F„(J)  will  still  be  a  Poisson  raindom  variable, 
but  with  parameter  where  refers  only  to  the  background  aerosol  component.  In  the 

ensuing  analysis,  it  is  assumed  that  {A(°^(j)}  do  not  depend  on  j,  viz.,  the  background  Poisson 
parameter  is  the  same  for  all  J  aerosol  samples  obtained.  Flirthermore,  it  is  assumed  that  is 

an  a  priori  known,  non-negative  sequence  due  to  the  presence  of  the  target  aerosol  signal. 

The  observed  count  <lata  {F„(J),j  =  1,2, ...,J}  is  a  set  of  mutually  independent  Poisson 
random  variables  with  parameters  if  the  signal  aerosol  is  present  (hypothesis  Hi) 

or,  with  parameters  A^®)(j)  =  A^°^  if  the  signal  aerosol  is  alxsent  (hypothesis  Ho).  This  forms  the  basis 
for  the  conventional  binary  hypothesis  testing  paradigm  in  which  the  system  receives  a  sequence  of 
aerosol  sample  count  data  F„(j)  from  J  independent  sampling  periods  that  can  correspond  to  a 
Poisson  process  with  parameter  A^;)  such  that 

Ho:X(j)  =  X<°\  j  =  l,2 . J; 

Hi  :  XU)  =  S(J)  +  X^°\  >  =  1,2 . J. 
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Detection  of  the  signal  aerosol  can  then  be  formulated  as  testing  the  null  hypothesis  Hq  versus 
the  alternative  hypothesis  Hi.  This  objective  can  be  achieved  with  the  likelihood  ratio  test  whose 
purpose  is  to  test  hypothesis  Hi  against  Ha,  holding  the  probability  that  Hi  is  accepted  when  Ha  is 
actually  true,  i.e.,  the  probability  of  false  alarm,  below  some  fixed  error  level. 


The  detector  is  based  on  a  generalized  likelihood  ratio  test  and  is  localized,  viz.,  it  examines 
a  set  of  count  data  spanning  a  relatively  short  sampling  duration.  While  not  structurally  necessary, 
localized  detection  is  better  suited  to  the  nature  of  the  signals  that  are  assumed  to  be  present  (time- 
varying  and  transient)  than  a  global  detection.  The  likelihood  ratio  A  is  the  ratio  of  the  probability 
(“likelihood”)  p({F„(/)}|/fi)  of  observing  the  count  data  {Fn(j‘),ji'  =  1,2,  ...,J}  under  Hi  to  the 
probability  of  observing  the  data  under  Hq,  p({^nO)}|^o)  so 

,  p{{F„(j)}\Hi) 

p{{FnU)}\Hoy 

If  is  known,  then  as  already  noted,  {Fn(j))  are  independent,  Poisson  random  variables,  so 
the  Joint  probability  p{{F„(j)}\Hi)  (t  =  0,1)  is  the  J-fold  pro<liict  of  the  individual  probability 
p(f’„(y)|//,)  of  obstirving  each  F„(J).  Consequently, 

p({F„(,)}|/fO  =  J]L-m_exp(_AWC;)) 


and  so  it  follows 


niSO)  +  A<»ll''-"'exp{-(SC))  +  Al«)) /f„y)l 


J=1 

Some  straightforward  algebraic  manipulations  can  be  performed  to  simplify  A  to  the  following  form: 


=  exp  (  -]^5C7) 

V  ]=i 


1  + 


fp(7) 


The  likelihood  ratio  must  now  be  compared  against  some  aj)propriatc  threshold  in  order 
to  determine  which  hypothesis  to  accept.  Since  A  is  used  by  the  d(x:ision  rule  only  in  comparison 
to  a  threshold — call  it  Aa — any  one-to-one  and  monotonic  transformation  on  A  may  be  [H:rformod 
on  both  A  and  Aa,  thereby  generating  a  modified  likelihood  ratio  to  be  compared  with  a  modifitxl 
threshold.  In  light  of  this,  it  is  convenient  to  consider  the  log-likelihood  ratio 


C  =  \nA 
.1 

=  E 

>=i 


FniJ) 


i„( 


(0) 
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In  its  present  form,  the  log-likclihood  ratio  detection  statistic  embodied  in  E^q.  (9)  has  been  derived 
under  the  implicit  assumption  that  is  known.  Since  it  is  unlikely  that  will  be  known  a 
priori,  we  must  now  confront  the  problem  of  estimating  A^®^  Intuitively,  it  should  be  expected  that 
if  the  target  aerosol  signal  is  absent,  then  A^®^  can  be  estimated  by  averaging  the  count  data  F„(j) 
over  the  sample  period  index  j.  Strictly  speaking,  in  the  event  that  A(°1  is  unknown,  the  likelihood 
ratio  test  requires  that  the  parameter  A^^l  be  obtained  as  the  maximum  likelihood  estimate  under  the 
null  hypothesis  Hq.  In  this  regard,  the  necessary  estimate  is  obtained  by  choosing  A^®^  to  maximize 
the  probability  of  receiving  the  count  data  samples,  viz.,  to  maximize 

PaFnU)mo)  =  n  -'>^1)!"^^  exp(-A(“)0))- 

The  maximum  likelihood  estimate  is  obtained  by  solving  the  equation  dp({F„(j)}\Ho) / dA^®)  =  0  and 
this  process  leads  to 


dA(«) 


=  0, 


from  which  it  follows 


Note  that  this  is  simply  the  average  of  the  count  data  F„(j)  for  the  case  when  the  signal  aerosol  is 
absent. 


Choosing  the  detection  threshold  is  just  as  important  as  computing  the  log-likelihood  ratio. 
It  is  important  to  note  that  the  choice  of  the  threshold  is  intimately  related  to  the  evaluation  of 
the  performance  of  the  detector,  and  to  this  purpose,  the  detection  threshold  T  must  be  chosen  to 
provide  a  given  false  alarm  rate  for  the  detector.  This  particular  choice  of  the  detection  threshold 
provides  the  maximum  probability  of  detection  for  a  given  probability  of  false  alarm.  It  is  important 
to  note  that  the  choice  of  the  “correct”  probability  ol  false  alarm  is,  of  course,  a  question  for  the  user 
of  the  detector  to  decide  and  one  on  which  physical  science  has  no  bearing.  However,  it  should  be 
stressed  that  it  is  not  possible  to  operate  any  detection  system  based  on  zero  false  alarm  probability; 
it  can  be  never  be  absolutely  guaranteed  that  a  zero  probability  of  false  alarm  will  occur  for  any 
physically  realizable  detector.  The  probability  of  false  alarm,  is  defined  as 

PvK  =  Pt{C  >  r|//o} 

and  the  associated  probability  of  detection,  Pp,  is  given  by 

Po  =  Pr{C>T\lli}. 
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The  calculation  of  the  detection  threshold  T  using  the  0cact  multivariate  Poisson  probabil¬ 
ity  distribution  is  a  diflicult  problem.  However,  for  sufliciently  long  sampling  times  Ut  the  count 
data  and  even,  more  pronouncedly,  the  detection  statistic  C  would  tend  to  Gaussian  vari¬ 

ates.  Hence,  a  Gaussian  approximation  for  C  allows  T  to  be  more  easily  computed.  Applying  the 
Lindelbcrg  condition  of  the  central  limit  theorem  [4]  to  the  sum  of  random  variables  in  Eq.  (9),  it 
is  straightforward  to  show  that  the  log-likelihood  ratio  can  be  approximated  by  a  Gaussian  vari¬ 
ate  whose  mean  and  variance  are  determined  as  follows.  In  the  absence  of  a  target  aerosol  signal 
(hypothesis  //q),  the  log-likelihood  ratio  C  has  a  mean  and  variance  given  by 

and 

Var(i:|;f„)=^ln»(l  +  ^)Al»>. 

>=1 

Similariiy,  in  the  presence  of  the  target  aerosol  signal  (hypothesis  //j),  the  mean  and  variance  of  the 
log-likelihood  ratio  assumes  the  form 

E(£|/fO  =  E(£l//o)  -h  ^  50)  In  (l  + 

and  ^ 

Var(£l//,)  =  Var(£|/fo)  +  ^  S{j)  In^  (l  +  ^)  • 

'J'hc  approximate  distribution  of  C  under  //q  enables  the  determination  of  the  tbreshold 
of  any  desired  probability  of  false  alarm.  The  approximate  distribution  of  C  under  Hy  enables  the 
determination  of  the  probability  of  detection  as  a  function  of  the  signal-to-noise  ratio.  Consequently, 
in  light  of  these  relationships  and  the  Gaussian  approximation  for  the  log-likelihood  ratio  C,  it  is 
evident  that  to  achieve  a  probability  of  false  alarm,  of  a  requires  the  solution  of 

Pfa  =  a  =  i  -  ierf(x/>/2),  (10) 

with  the  detection  threshold  calculated  from 

T  =  E[£|//o]  -»-x^Var[£|//o].  (11) 

Here,  erf(')  denotes  the  error  function.  For  this  choice  of  detection  threshold,  the  probability  of 
(h!tection  can  then  be  compiiteel  from 

/’u  =  ^  -  ^<’rf(y/\/2),  (12) 
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where 

y={T-  E[/:|tf,])  j ^Var[i:|/f,]. 


(13) 


For  the  likelihood  ratio  test  that  has  been  formulated,  the  performance  of  the  detection 
statistic  depends  only  on  the  dimensional  parameter  J,  the  Poisson  parameter  associated  with 
the  background  aerosol  samples,  and  the  signal-to-noise  ratio  (SNR)  which  is  defined  as 


In  examining  the  probability  of  detection,  there  are  several  different  parameters  that  can  be  varied. 
In  the  following  simulations,  a  Poisson  parameter  A^®)  =  10  is  used  for  the  representation  of  the 
aerosol  background  noise  level.  In  Figures  7-8,  the  probability  of  detection  (fb)  is  shown  against 
the  probability  of  false  alarm  (/Va)«  parametrized  by  the  SNR  for  J  =  10.  Recall  that  J  is  the 
number  of  aerosol  samples  used  in  the  computation  of  the  log-likelihood  ratio.  Figures  9-10  are 
analogous,  but  they  correspond  to  the  case  for  J  =  25  aerosol  samples  used  in  the  detection  process. 
As  would  be  expected,  observe  that  the  greater  J  is,  the  better  the  performance  of  the  detection 
algorithm.  These  curves  show  the  fact  that  for  a  fixed  SNR,  the  detector  with  the  most  aerosol 
samples  has  the  best  detectability.  Figures  11  and  12  depict  the  probability  of  detection  versus 
t  he  signal-to-noise  ratio  for  several  prescribed  (i.e.,  fixed)  probabilities  of  false  alarm  for  J  =  10 
and  25,  respectively.  It  is  interesting  to  note  that  for  SNRs  greater  than  or  equal  to  about  0.4, 
the  detection  performance  curves  suggest  that  PpA  depends  only  on  the  number  of  aerosol  samples, 
J ,  used  in  the  detection  process  and  is  independent  of  the  SNR  (viz.,  of  the  signal  level  of  the 
t^u■get  aerosol) .  Consequently,  in  this  regime  of  operation,  the  detection  threshold  can  be  calculated 
without  knowledge  of  the  signal-to-noise  ratio  and  need  not  be  adjusted  to  keep  a  constant  false 
alarm  rate  throughout  the  detection  procedure.  In  this  regime,  the  detection  algorithm  exhibits 
constant  false  alarm  rate  (CFAR)  behavior. 


In  summary  then,  the  log-likelihood  ratio  detection  algorithm  can  be  implemented  as  follows: 
(1)  choose  a  for  the  desired  Pvk  using  Eqs.  (10)  and  (11);  (2)  compute  the  background  Poisson 
parameter  A^°^  using  only  samples  obtained  from  the  pure  aerosol  background;  (3)  compute  the 
log-likelihood  ratio  £  (cf.  Eq.  (9));  (4)  calculate  the  detection  threshold 't  from  Eqs.  (12)  and  (13); 
and,  (5)  compare  £  with  T  and  declare  the  target  (signal  aerosol)  to  be  present  if  £  >  T. 
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IV.  SOME  NUMERICAL  EXAMPLES 

In  this  section,  the  use  of  the  detection  procedure  based  on  the  log-likelihood  ratio  statbtic, 
is  illustrated  by  some  numerical  examples.  The  background  aerodynamic  particle  size  distribution 
function  is  assumed  to  be  represented  by  a  gamma  size  distribution  with  the  shape  parameter  a  =  1 .5 
and  the  scale  parameter  fi  =  2.5  /tm.  In  the  following  examples,  the  total  number  of  aerosol  particles 
in  the  aerodynamic  diameter  range  [Di,D2]  =  [1.0,15.0]  /trn  in  the  sampltxl  volume  i.s  n.s*Ml  as 
the  count  data  for  the  construction  of  the  log-likelihood  ratio  detection  statistic. 

Example  1:  Figure  1.3  shows  a  target  aerosol  signal  in  the  form  of  the  time  history  of  the 
number  concentration  at  a  receptor  position  where  an  APS  has  been  deployed.  In  this  example, 
the  peak  amplitude  of  the  number  concentration  is  40  cm“® — the  time  of  passage  of  the  quasi- 
instantaneous  puflF  over  the  receptor  position  is  about  15  min.  For  this  example,  the  particle  number 
density  N*  of  the  gamma  size  distribution,  used  to  model  the  background  particle  size  distribution, 
was  adjusted  to  yield  a  background  aerosol  count  with  an  intensity  =  10.0.  A  sampling  time 
of  4  =  60  s  was  chosen  for  the  APS  detection  system  operating  at  the  aerosol  flowTate  Q,  =  8.3..3 
cm*  s~’ — this  sampling  time  (i.e.,  exposure  time)  was  chosen  to  ensure  a  statbtical  sampling  error 
of  smaller  than  or  equal  to  0.5  percent.  The  detection  performance  will  be  expressed  in  terms  of 
the  probability  of  detection,  Pq,  as  a  function  of  the  number  of  acro.sol  samples,  ./,  use<l  in  the 
detection  scheme,  for  a  fixed  probability  of  false  alarm  Pka-  For  this  example,  Pyj^  =  10“^,  10"*,  aiul 
10~*  were  chosen.  It  is  assumed  that  aerosol  samples  used  in  the  detection  procedure  were  collected 
immediately  after  the  initial  arrival  of  the  transient  target  aerosol  signal.  The  performance  of  the 
detector  b  shown  in  Figure  14.  Observe  that  the  probability  of  detection  Pd  m  1.0  for  J  >  7  and 
that  the  detection  probability  is  independent  of  Pp^  after  this  point.  With  the  prescribed  sampling 
time  of  ^  =  60  s,  this  implies  that  the  APS  was  able  to  detect  the  target  signal  with  certainty 
approximately  7  minutes  after  the  initial  arrival  of  the  transient. 

Example  2;  The  experiment  described  in  Example  1  was  repeated,  but  this  time  the  target 
aerosol  signal  piossesses  the  form  depicted  in  Figure  15.  The  target  signal  is  a  broader,  semi- 
continuous  puff  with  a  peak  amplitude  of  approximately  1.4  cm"* — the  time  of  passage  of  this 
transient  signal  over  the  receptor  position  is  about  150  min.  For  this  example,  the  particle  number 
density  N'  was  adjusted  to  provide  a  background  intensity  of  =  5.0.  It  is  noted  that  the 
background  aerosol  clutter  level  is  lower  than  that  as.sumed  in  F/Xampio  1.  A  sampling  time  t,  =  20 
s  was  chosen  for  the  APS,  opt'ral.iiig  at  the  volume  flowrate  of  Q,  =  S.3..3  cm*  s  '  in  onler  to 
yield  a  prescribed  value  for  the  relative  fluctuation  intensity  of  the  total  number  of  aerosol  partich's 
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in  the  sampled  volume  of  1  percent  (viz.,  this  means  that  one  can  expect  a  standard  deviation 
in  the  total  number  of  particles  of  approximately  I  percent  of  the  expected  total  number  due  to 
Poisson  distributed  fluctuations  of  the  number  of  particles  sampled  about  their  mean  gamma  size 
distribution).  For  this  example,  the  number  of  samples  used  for  the  detection  procedure  is  /  =  45 
and,  with  a  sampling  time  of  ^  =  20  s,  this  implies  that  a  detection  decision  is  made  after  every 
sampling  duration  of  15  min.  Figure  16  shows  that  detection  probability,  at  prescribed  fabe  alarm 
probabilities  of  10~®,  10“^  and  10“^,  for  each  sampling  duration  of  15  min  after  the  initial  arrival  of 
the  target  signal.  In  the  plot,  the  time  index  along  the  absebsa  corresponds  to  the  number  index 
of  each  detection  event  of  sampling  duration  15  min.  In  this  case,  the  detection  performance  b  fair 
since  for  /Va  =  0.01,  there  is  a  0.87  probability  of  detection  of  the  target  signal  (shown  in  Figure  15) 
for  the  detection  period  that  includes  the  peak  of  the  target  signal. 


V.  CONCLUSIONS 

In  thb  paper,  the  basic  theory  and  algorithms  for  a  statistical  sampling  error  model  and  a 
detection  scheme  for  the  aerodynamic  particle  size  analyzer,  has  been  developed.  It  has  been  shown 
how  the  Poisson  process  can  be  used  as  the  basb  for  a  sampling  error  model  which  characterizes 
the  fluctuations  of  a  particle  size  dbtribution-related  variable  due  to  stochastic  fluctuations  of  the 
number  of  particles  sampled  in  each  particle  size  interval.  The  model  developed  here  allows  for  the 
proper  specifleation  of  the  required  sampling  program  for  the  APS  in  order  that  the  sampling  error 
he  below  some  flxed  error  level.  With  the  prescription  of  an  appropriate  sampling  program,  this 
paper  presents  an  algorithm  for  detecting  a  target  aerosol  signal  against  a  clutter  aerosol  background. 
I'he  detection  algorithm  is  based  on  the  gcticralize<l  likelihood  ratio  test  and  can  be  applied  to  the 
detection  of  an  aerosol  signal  sequence  in  a  set  of  J  aerosol  samples  with  a  common  aerosol  clutter 
background.  It  has  been  shown  how  the  probability  of  detection  of  the  detector  can  be  determined 
for  a  prescribed  probability  of  false  alarm,  including  the  effects  of  the  signal-to-noise  ratio,  the 
intensity  of  the  background  aerosol  clutter,  and  the  number  of  aerosol  samples,  J,  used  in  the 
detection  process.  Some  detection  results  have  been  presented  which  illustrate  the  wide  variations 
in  performance  which  may  be  observed,  dependent  on  the  background  clutter  condition,  the  form  of 
the  target  aerosol  signal  to  be  detected,  the  choice  of  the  detection  threshold,  and  the  signal-to-noise 
ratio. 


Practical  application  of  the  detection  algorithm  described  in  this  paper  must  also  consider 
the  issues  related  to  the  operational  characteristics  of  the  APS.  One  issue  concerns  the  efficiency 
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of  the  APS,  viz.,  what  fraction  of  the  particles  in  the  sampled  volume  are  actually  counted  by 
the  instrument.  Another  issue  is  the  robustness  of  the  APS  and  the  concomitant  failure  rate  of 
the  instrument  when  employed  in  a  continuous  held  sampling  mode  for  monitoring  and  detection. 
Finally,  the  paper  has  only  addressed  the  detection  performance  of  a  single  APS  operating  as  a 
centralized  system — the  design  and  performance  of  a  distribuUKi  signal  detection  system  based  on 
an  array  of  aerodynamic  particle  size  analyzers  which  incorporates  the  possibility  of  sensor  failure 
needs  to  be  addressed. 
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FIGURE  4 

Curve  of  (N*  VJ)”''*/  as  a  function  of  the  gamma  size  distribution  shape  parameter  a  for  ^  =  1.0  /tm  and 
n  =  3. 
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FIGURE  8 


Probability  of  detection  (/b)  versus  probability  of  false  alarm  (FV-a)>  shown  for  several  signal-to-noisc  ratio 
(SNR)  values,  given  J  =  10. 
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FIGURE  9 

Probability  of  detection  (/^)  versus  probability  of  false  alarm  (/Va),  shown  for  several  signal-to-noisc  ratio 
(SNR)  values,  given  J  ~  25. 
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FIGURE  10 

Probability  of  detection  (fb)  versus  probability  of  false  alarm  (Pfa),  shown  for  several  signal-to-noise  rat  io 
(SNR)  values,  given  J  =  25. 
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FIGURE  11 

I’robability  of  dctocUon  {I’d)  versus  signal-to-noise  ratio  (SNR)  for  Ppj^  =  10  '®,  10'  10  10  ’’  and  ./ 
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FIGURE  14 

Probability  of  detection  (Pp)  for  the  transient  aerosol  signal  depicted  in  Fig.  13  as  a  function  of  the  number 
of  aerosol  samples,  J,  used  in  the  detection  process. 
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FIGURE  15 

Time  history  of  the  number  concentration  of  the  target  aerosol  signal  used  in  Example  2. 
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FIGURE  16 

Probability  of  detection  {Pq)  for  the  transient  aerosol  signal  depicted  in  Fig.  15  as  a  function  of  the  detection 
period  (indicated  by  the  time  index).  Each  detection  period  consisted  of  a  sampling  duration  of  15  min. 
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The  deteniination  of  the  distribution  of  airborne  toxic  particles  as  a  function  of 
the  aerodynaaic  diaaeter  provides  iaportant  inforaation  as  well  as  criteria  for  the  def¬ 
inition  of  haxard  as  applied  to  levels  of  airborne  contaainatftjn.  This  is  because  the 
aerodynaaic  particle  size  distribution  eabodies  the  inforaation  related  to  particle 
density,  diaaeter,  shape  factor  and  slip  correction  that  is  critical  for  the  character¬ 
ization  of  particle  mo^ftion  in  settling  and  iapaction  and  it  is  these  actions  that  are 
responsible  for  particle  deposition  in  the  respiratory  tract  and  particle  collection  in 
aerosol  saapling  devices.  For  a  given  definition  of  hazard  based  on  some  paraaeter 
related  to  the  aerodynaaic  size  distribution,  this  paper  develops  a  statistical  saapling 
error  aodel  for  the  paraaeter  that  is  based  on  the  Poisson  process.  Given  that  an 
appropriate  sampling  program  has  been  designed  for  the  measurement  of  the  size  distri¬ 
bution-related  paraaeter  with  the  aerodynaaic  particle  size  analyzer,  this  paper  proceed] 
to  the  derivation  of  an  optisiua  detection  algorithm  for  the  detection  of  a  signal 
aerosol  sequence  in  a  set  of  J  aerosol  samples  with  a  coaaon  background.  The  detection 
algorithm  is  based  on  the  generalized  likelihood  ratio  test  in  which  the  received  count 
associated  with  the  aerosol  sample  is  modeled  as  a  Poisson  distributed  random  variable. 
Performance  analyses  of  the  resulting  algoritlui,  based  on  the  probability  of  detection 
(P  )  versus  signal-to-noise  ratio  for  several  given  fixed  false  alarm  probabilities 

(P  ),  are  presented. 
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